Parametric Amplification of Gravitational Fluctuations During Reheating 
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We demonstrate that cosmological perturbations can undergo amplification by parametric reso- 
nance during the preheating period following inflation, even on scales larger than the Hubble radius, 
without violating causality. A unified description of gravitational and matter fluctuations is crucial 
in order to determine the strength of the instability. To extract specific signatures of the oscillating 
inflaton field during reheating, it is essential to focus on a variable describing metric fluctuations 
which is constant in the standard analyses of inflation. For a massive inflaton without self coupling, 
we find no additional growth of super-horizon modes during reheating beyond what the usual anal- 
yses of the growth of fluctuations predict. However, for a massless self coupled inflaton, there is an 
initial increase in the value of the gravitational potential, which may lead to different predictions 
for observations. 
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1. As was initially realized in Q and worked out in 
detail in |^]-^| and many other papers |^] , parametric res- 
onance instabilities play a crucial role in the dynamics 
of the reheating of an inflationary Universe. Reheating 
(or now more accurately called "preheating" 0) is the 
period after inflation when the inflaton, the scalar field 
driving the period of exponential expansion, oscillates 
coherently about its ground state and gradually trans- 
forms its energy into the matter and radiation which the 
present Universe is made up of. 

The parametric resonance instability can be seen by 
considering the linearized equation of motion of fields x 
which couple to the inflaton. Neglecting for a moment 
the expansion of the Universe, the equation of motion 
for the Fourier modes Xk becomes a harmonic oscillator 
equation with a periodically varying mass, the Mathieu 
equation. It is well known that this equation admits in- 
stability bands, regions of k for which the solutions grow 
exponentially. As already noted in jl] || and discussed 
in detail in the instabilities persist after taking the 
expansion of the Universe into account. 

The parametric resonance instabilities have important 
consequences for cosmology. They will lead to a reheat- 
ing temperature which can be much larger than would be 
obtained by calculating the efficiency of reheating using 
perturbation theory, as was done initially in j^]]. This 
could have important implications for GUT-scale baryo- 
genesis Q . Since parametric resonance only amplifies the 
amplitude of certain frequency bands of excitations, the 
state produced will initially have a nonthermal frequency 
distribution. This may lead to the possibility of the for- 
mation of topological defects after inflation j| |l3 , and 
to the production of supermassive dark matter |L5fl . 

Initially (l), the parametric instability was discussed 
in a toy model of new inflation in which the resonance 
bands were narrow. It was then pointed out ^ that in 
more realistic models of chaotic inflation, the instability 
bands were in fact much wider ( "broad resonance" ) , and 



the Floquet exponent fik giving the rate of exponential 
growth correspondingly larger. In |l6| it was discovered 
that a negative sign of the coupling constant between the 
inflaton and the x field leads to a greatly enhanced insta- 
bility ( "negative coupling instability" ). White noise elim- 
inates the stability bands alltogether |l7];[l8| and mode by 
mode strengthens the resonance. 

The inflaton field also couples to linearized metric per- 
turbations. Bassett [|l9| first realized that gravitational 
waves also experience parametric resonant amplification 
during reheating. Scalar metric perturbations, however, 
constitute the components of the metric which couple di- 
rectly to the matter sector. Hence, it is to be expected 
that scalar perturbations are amplified more strongly 
than gravitational waves. The growth of scalar perturba- 
tions due to the oscillations of the inflaton field was first 
considered by Nambu and Taruya pp|] and by Kodama 
and Hamazaki |2l]] (see also ^2[). Nambu and Taruya 
concluded that scalar perturbations are amplified during 
reheating, but did not compare their growth with the 
usual growth of cosmological perturbations. Kodama and 
Hamazaki focused on the evolution of the "Bardeen pa- 
rameter" , a gauge invariant measure of the cosmological 
perturbations which in the usual analysis of the growth 
of fluctuations (which neglects the oscillations of the in- 
flaton field) is constant in time for modes with wave- 
length larger than the Hubble radius. They concluded 
that in spite of the singular behavior of the the quantity 
c s = Pi P- where p and p are pressure and energy density, 
respectively, the Bardeen parameter remains constant. 
On the other hand, in a very interesting recent letter, 
Bassett, Kaiser and Maartens ||^] have re-analyzed this 
problem and argue that there is a negative coupling para- 
metric resonance instability which leads to a rapid growth 
of metric perturbations which in turn act as a pump for 
matter perturbations. Consequently, matter perturba- 
tions are amplified parametrically even in matter models 
in which there is no resonance in the absence of gravita- 
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tional fluctuations. 

In this letter, we analyze the growth of metric inho- 
mogeneities during reheating in a more complete way, 
making use of the gauge-invariant theory of perturba- 
tions (see p4| for a review). Since matter and metric are 
coupled by the Einstein constraint equations, the fluctu- 
ations can be described completely by a single gauge- 
invariant variable <£>. In longitudinal gauge, the per- 
turbed metric can be written in terms of $ as 



ds 2 = dt 2 (l + 2$) - a 2 (t)(l - 2$) 



(1) 



where a(t) is the scale factor. The equation of motion for 
<f> (linearized Einstein-Higgs equation) is a linear second 
order differential equation. As pointed out in (2^] , two of 
the coefficients in this equation of motion are singular due 
to the oscillations of the inflaton field. This is a problem 
encountered a long time ago jE5[in the context of axion 
fluctuations. As realized in [p0||21| , the divergence disap- 
pears if, instead of <E>, one considers the equation of mo- 
tion for the Mukhanov variable Q |p6| , a variable in terms 
of which the quantization of cosmological perturbations 
is straightforward (see for a review). We demonstrate 
that an instability persists in the equation of motion for 
a rescaled variable Q. This instability, however, is not of 
negative coupling type. For a massive inflaton, it only 
leads to an increase of Q proportional to a 3//2 (t) for long 
wavelength fluctuations. Hence, the amplitude of Q is 
constant in time, and there is no amplification of fluc- 
tuations beyond what the usual theory (which does not 
take the details of reheating into account) predicts. How- 
ever, for a massless self coupled inflaton, Q experiences 
an initial increase. From the regular behavior of Q, a 
completely regular behavior of $ can be deduced. We 
show why the usual methods to study the evolution of 
perturbations in inflationary cosmology miss additional 
growth of fluctuations due to the oscillating inflaton field. 

2. Our starting point are the equations of motion for 
the perturbations of the Einstein-Higgs system about a 
Friedmann-Robertson- Walker background solution. In 
terms of the gauge invariant metric and matter variables 
4> (see (jl|)) and 6<j) g i (which in longitudinal gauge is equal 
to the scalar field perturbation 5<fi, because of which we 
will subsequently drop the subscript), the system of equa- 
tions in momentum space is 

k 2 

$ + 3ff$ + [— + 2(H + H 2 )} $ = K 2 {<t> + H<j))8(t> (2) 



3H5<f> + (-j + V")5(j) = 4$0 - 2V'$ (3) 



$ + = -k 2 <M</>, (4) 



field with potential energy density V ((/)), and where a 
prime denotes the derivative with respect to <p. Equation 
(3) is the equation of motion for S(f> in a perturbed metric, 

(1) is the Einstein momentum constraint equation, and 

(2) is a combination of the dynamical equation of motion 
for $ and the Einstein energy constraint equation (see 
eqs. (6.42) and (6.40) of Q). These are exactly the 
same equations (except expressed in terms of physical 
time) as eqs. (2 - 4) of |||. 

Because of the Einstein constraint equation ([|), there 
is only one physical degree of freedom which we can 
choose to be 4>. Note that the source term in (||) is not 
suppressed compared to the terms on the left hand side of 
the equation in spite of the factor of k 2 which multiplies 
the term, since the constraint equation involves a com- 
pensating factor of k 2 . The correct equation of motion 
for $ is obtained by inserting the constraint equation (^) 
into (0), with the result 



3> + (H - 2-^)$ + (— + 2H - 2HX-)<S> = . 



(5) 



During the slow rolling period of an inflationary cos- 
mology, the coefficients in this equation are well-behaved. 
However, oscillations of <fi during reheating lead to sin- 
gularities. As was realized in p0| , pl| , these singularities 
can be eliminated by making use of Mukhanov's variable 
p6[ Q, the combination 



Q 



H 



(0) 



of the gauge invariant matter and metric perturbations in 
terms of which the unified quantization of the matter and 
metric perturbations is easy. In terms of Q, the equation 
of motion (|) becomes [||j20||2l]] 



Q + 3HQ + (V" + % + 2{^-+ 3H))Q = . (7) 

As is evident, the coefficients of this differential equation 
are regular. Given Q, it is possible to obtain $ since (|^) 
can be rewritten in the form 



a 2 2 H [ 6^ ) 



The Hubble damping term in the equation (Q) for Q 
can be eliminated by introducing the rescaled variable Q. 
For a massive inflaton with potential V(4>) — m 2 <j) 2 /2, 
Q = a 3 / 2 Q. In terms of Q, (@) becomes 



Q + [V" 



2( 



H 
H 



w)-\(h 2 + \h)]q = q. 

(9) 



where k 2 = 8irG, H = a/a is the Hubble expansion rate, 
4> is the homogeneous background field for a scalar matter 



Making use of the background Einstein equations, (g) can 
be written as 
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2k 
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H 



+ 3k 
2H 2< 
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(10) 



3k 2 

+ —rPAQ = °' 



where is the background pressure of the scalar field. 

After the period of slow rolling has ended, the value 
of H is smaller than m. Hence, it follows from the back- 
ground equation of motion for (f> that - in the absence of 
back-reaction and with accuracy increasing in time - the 
motion of 4> = a 3 / 2 cf> is oscillatory in time. In this limit, 
equation (pf) has the form 



where 



Q + [A(k) - 2gcos(mt)] Q = , 



A(k) = m 2 + 



(11) 



(12) 



where r contains the time average of the last four terms 
in the square bracket of (|Io|), and q contains the coeffi- 
cients of the oscillating parts of these terms. Since q(t) 
is decreasing, ( |li"|) is not of the form of the usual para- 
metric resonance equation, and no exponentially growing 
solutions will result. 

The second of the four last terms in the square bracket 
of ( [To| ) is the most important. Its initial amplitude is the 
largest, and it decays the least fast as a function of time. 
Approximating a{t) ~ t 2 / 3 corresponding to a pressure- 
less phase, it is easy to check that the second term decays 
as t~ l , whereas the other three terms decay as t~ 2 . Note 
that the decay rate of the dominant term of q as a func- 
tion of time is less fast than the corresponding decay rate 
of q for matter fluctuations 0], a point already empha- 
sized by p3[ . The amplitude of this leading term in q 
starts out slightly larger than m 2 . Due to the expansion 
of the Universe, the instability does not lead to exponen- 
tial increase in Q, but only to an increase proportional to 
a 3 / 2 (t) (see Figure 1), which implies that the amplitude 
of Q remains constant. 

The results are different for a massless self coupled 
inflaton with potential V((f>) — A0 4 /4. In this case, we 
eliminate (following f27|l ) the Hubble damping term in 
(0) by introducing the rescaled variable Q = a{rf)Q and 
by working in terms of conformal time rj. In this case, 
(0) becomes 

Q+ [A(fc) + 3c7i 2 (x,-U]g =0, (13) 
v2 

where x = s/XA^T] is a rescaled conformal time, A$ is 
the amplitude of oscillation of a(rj)(p, cn stands for the 
elliptic cosine function, and 



A(k) 



(14) 




Here, s stands for oscillatory terms which are important 
initially and lead to an increase in Q (even for k 2 <C 
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FIG. 1. Evolution of Q as a function of time in the case 
of a massive inflaton for a mode with k — 0.1m, with initial 
conditions such that Q is — 1 and its time derivative at the 
initial time when <f> = 0.2m p i, after slow rolling has ended. 
Time is expressed in units of m _1 . 



A.4H W ), but which decay in time as 77 1 . For s = 0, 
Equation (13|) is a Lame equation, has been studied in 
detail in p7[ and exhibits no resonance for k 2 -C XA 2 hi . 

3. We have solved the equation of motion ([[o|) numer- 
ically. Figure 1 shows the resulting time evolution of the 
rescaled variable Q in the case of a massive inflaton for 
the mode k = 0.1m (about six times larger wavelength 
than the Hubble radius) over a period of several oscilla- 
tions of the background field. As mentioned above, the 
amplitude of Q grows linearly in time (and thus the am- 
plitude of Q remains constant) to within the numerical 
accuracy. 

Does the growth of Q imply that there is an increase in 
the amplitude of cosmological fluctuations beyond what 
is calculated in the usual analyses which model the re- 
heating as a continuous change in the equation of state 
without any oscillations? According to the usual treat- 
ment (see e.g. Q), the details of the equation of state 
are completely unimportant for the final amplitude of 
fluctuations on scales larger than the Hubble radius dur- 
ing the reheating period. If there were an increase in 
the amplitude of $ beyond what is usually predicted, 
this could have important implications for the spectrum 
of density perturbations and microwave anisotropies, as 
conjectured in p3| . First of all, let us note that such an 
effect would NOT violate causality. Inflation has already 
set up fluctuations on all scales smaller than the Hubble 
radius at the beginning of inflation, and in modelling the 
Universe as a Friedmann Universe, we have inserted cor- 
relations on even larger scales. These perturbations can 
self gravitate and increase in magnitude during reheating 
without violating causality. This is precisely the way in 
which quantum fluctuations created during inflation be- 
come large after inflation. Since the phase of oscillation 
of 4> is coherent over a region much larger than the Hub- 
ble radius during reheating, it is therefore possible that 
the oscillations might have an effect on modes larger than 
the Hubble radius. 

The clearest way to determine if the growth of Q cal- 
culated above is a new effect is to calculate the quantity 
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(15) 



which according to the standard theory of the growth of 
cosmological perturbations remains constant for modes 
with wavelength larger than the Hubble radius j28j . Here, 
w = pj p describes the equation of state, and p and p 
are the pressure and energy density, respectively. More 
precisely j24j], (j> 2 (, = is equivalent to the equation of 
motion for $ (for modes with wavelength larger than the 
Hubble radius). Hence, it is usually deduced that £ = 0. 
This conclusion, however, may break down if (f> = which 
occurs precisely in the phase we are studying here dur- 
ing which <j> is oscillating coherently. Hence, it is possi- 
ble that additional resonant amplification of fluctuations 
during reheating occurs. 

The variable £ is related to Q via 



H 



(16) 



Since the amplitudes of H and (f> decrease at the same 
rate (proportional to t^ 1 ) if we assume that a(t) ~ i 2 / 3 , 
the constancy of the amplitude of Q (for a massive in- 
flaton) leads to the conclusion that £ is constant when 
evaluated at the same phase during each oscillation pe- 
riod of the inflaton <fi. Since l + w has constant amplitude 
during the period of oscillation, it follows from the con- 
stancy of £ that the amplitude of <£> is constant, inasfar 
as the <3? term in ([l5[) can be neglected. This can also 
be seen by evaluating $ directly using (||) (see ^9|) for 
further discussion). 

Note that $ is the basic physical quantity which is 
well-defined at all times. It is the quantity which deter- 
mined the power spectrum of density fluctuations and of 
CMB anisotropics. In contrast, £ is an auxiliary quan- 
tity. At each zero of (f> there is a singularity in the re- 
lation between $ and £. What is therefore important 
is to calculate the value of £ for each zero crossing of (f> 
and compare the values. In Figure 2 we plot (1 + w)£ as 
a function of time, determined directly from (|l6|). The 
fact that the amplitude of oscillation of this function is 
constant implies that £ does not change over a period. 
This demonstrates that the growth of Q observed in Fig- 
ure 1 is in exact agreement with the usual analysis of the 
growth of cosmological perturbations. 

4. To summarize, we have studied the growth of cos- 
mological perturbations during a phase of coherent oscil- 
lations of the scalar field driving inflation, which leads to 
singularities in c 2 s = pj p. In this period, the equations of 
motion of the fluctuations are similar in structure to the 
equations describing coupled matter fields during reheat- 
ing. As emphasized in P, |l6| , p7t , the presence or absence 
of a mass in the inflaton field is crucial. To study the evo- 
lution of cosmological perturbations, it is crucial to work 
in terms of variables in terms of which the singularities 
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FIG. 2. Evolution of (1 + w)C, for a theory with a massive 
inflaton, and for same mode, initial conditions and initial time 
as in Fig. 1. 



are absent, and in which the usual growth of cosmolog- 
ical fluctuations (obtained without taking into account 
the oscillations of <fr) is factored out, such as £. 

For a massive inflaton, we find no amplification of long 
wavelength fluctuations (wavelength larger than the Hub- 
ble radius during reheating) - the amplitude of £ is con- 
stant. A modelling of the reheating period including the 
oscillations of <fi will lead to the same growth as is ob- 
tained in the usual analyses of perturbations in which 
the transition from the inflationary phase to the post- 
inflationary radiation-dominated period is modelled (im- 
plicitly) by a monotonous transition in the equation of 
state. In a companion letter |3C(| , Parry and Easther 
demonstrate that even a full nonlinear analysis does not 
lead to any additional growth of fluctuations with wave- 
length larger than the Hubble radius during reheating. 

However, for a massless self coupled inflaton, initially 
both the gravitational potential $ and £ grow for modes 
larger than the Hubble radius during reheating, as con- 
jectured in (see [^9| for a detailed discussion). 

Note that the negative coupling instability discussed 
in is not present in the relevant equation (Q) for Q. 
The negative coupling instability in the equation for <f> 
is precisely the instability which is responsible for the 
amplification of the fluctuations in the standard analysis 
of the growth of cosmological perturbations. 
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